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Abstract. The local electronic structure of YBa2Cu307 has been calculated using first-principles cluster 
methods. Several clusters embedded in an appropriate background potential have been investigated. The 
electric field gradients at the copper and oxygen sites are determined and compared to previous theoretical 
calculations and experiments. Spin polarized calculations with different spin multiplicities have enabled a 
detailed study of the spin density distribution to be made and a simultaneous determination of magnetic 
hyperfine coupling parameters. The contributions from on-site and transferred hyperfine fields have been 
disentangled with the conclusion that the transferred spin densities essentially are due to nearest neighbour 
copper ions only with marginal influence of ions further away. This implies that the variant temperature 
dependencies of the planar copper and oxygen NMR spin-lattice relaxation rates are only compatible with 
commensurate antiferromagnetic correlations. The theoretical hyperfine parameters are compared with 
those derived from experimental data. 
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1 Introduction 

Although discovered M in 1986, the high temperature su- 
perconducting (HTS) copper oxides continue to present 
challenge to both experiment and theory. The underlying 
mechanism (or mechanisms) responsible for their super- 
conductivity is still ambiguous. To understand the vari- 
ous contributions of magnetism and phonons in high tem- 
perature superconductivity requires a detailed knowledge 
of the electronic structure of HTS materials. These ex- 
hibit considerable structural complexity but this complex- 
ity brings the advantage of various nuclear species located 
throughout the unit cell which serve as local probes of 
the electronic and magnetic properties. Through nuclear 
quadrupolc and magnetic resonance (NQR and NMR) 
techniques, these nuclei probe passively the local elec- 
tronic structure and have delivered a tremendous amount 
of data about both static and dynamic properties of HTS 
materials M. 

The nuclei experience the spin polarization of the 
electronic system through various hyperfine interactions 
which basically are well known. It turned out, however, 
that in HTS copper oxides their analysis is not as straight- 
forward as expected. For example, the spin shifts which 
give information on the static spin susceptibility have to 
be separated from the chemical shifts which, being inde- 
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pendent on temperature, are assumed to dominate the 
Knight shifts at low temperature. Precise measurements 
of these NMR shifts in the superconducting state, how- 
ever, are difficult due to the field inhomogeneities pro- 
duced by the vortex lattice. Furthermore, the spin shifts 
at the planar copper nuclei are extremely small in all 
cuprates when measured with the applied field perpen- 
dicular to the planes. 

Dynamic properties of the spin fluid can be probed by 
the various nuclear relaxation rates. The interaction en- 
ergy is determined by the sum of the hyperfine interactions 
of which the transferred hyperfine fields make a crucial 
contribution as has first been pointed out in Refs. P,B- 
In particular, the variant temperature dependence oithc 
spin-lattice relaxation rates at planar copper and oxygen 
nuclei could be explained by assuming that the trans- 
ferred fields at 17 from the two neighbouring copper ions 
cancel if these moments are antiferromagnetically coupled 
whereas the copper nuclei could reflect the antiferromag- 
netic spin fluctuations. Neutron scattering measurements, 
however, indicated that the fluctuation peaks were incom- 
mensurate and not located at the reciprocal wave vectors 
Q = (7r/o,7r/o). 

Zha, Barzykin, and Pines H then advocated trans- 
ferred hyperfine fields at the oxygens from next nearest 
neighbour copper ions which would further suppress the 
contributions from the two nearest neighbours. With this 
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extended model they could explain a variety of NMR data 
in good agreement with incommensurate antiferromag- 
netic fluctuations as exhibited by neutron scattering re- 
sults. 

Although the focus of NQR and NMR investigations of 
HTS materials has recently changed to inhomogeneities, 
in particular static and dynamic stripe phenomena, some 
of the "older" problems still wait for definitive answers. 
Among these problems the applicability of a one-band ap- 
proach to describe the low-lying spin and charge excita- 
tions in these materials is of special importance. Several 
experiments PJ7fl9] indicate that a second spin fluid with 
independent dynamics is required. 

Despite the rich information these experiments pro- 
vide, there exist only few theoretical first-principles ap- 
proaches which address the determination of electric field 
gradients (EFGs) and magnetic hyperfine interactions. For 
the YBa2Cu306+a: system, EFGs at various nuclear sites 
have been early obtained by Das and co-workers p ]l(|[ll| 
with ab initio cluster calculations using the unrestricted 
Hartree-Fock (UHF) method-and by Schwarz and co- 
workers |]3|l3| and Yu et al. Q , who employed the full- 
potential linear augmcnted-planc-wave method (FLAPW) 
within the local density approximation (LDA). Results of 
Hartree-Fock (HF) calculations have also been published 
by Winter p5) . The EFGs calculated with these three dif- 
ferent methods agree, more or less, and reproduce apart 
from one exception the experimental data quite satisfac- 
torily. The exceptional case is the EFG at the planar Cu 
sites. Using a large cluster that contained 74 atoms, we 
performed |1(| calculations with the density functional 
(DF) method and obtained values for the Cu EFGs which 
were in better agreement with the experiments. These 
calculations were spin non-polarized and no information 
on magnetic hyperfine interactions was available. Spin 
polarized calculations, however, have recently been per- 
formed |17J for the La2Cu04 system which allowed the 
study of hyperfine coupling parameters. 

In this paper, we report the results of extensive cluster 
studies of the electronic structure of YB^CusOy. Spin po- 
larized calculations with the DF method with generalized 
gradient corrections to the correlation functionals have 
been performed for various clusters. The resulting elec- 
tronic structure, the charge and spin distributions have 
been analyzed. In particular, the on-site and transferred 
hyperfine fields have been investigated in detail and it is 
demonstrated, that transferred interactions from copper 
ions beyond the nearest neighbours are small for both pla- 
nar copper and oxygen nuclei. The relevance of this finding 
for the interpretation of 17 NMR spin-lattice relaxation 
is discussed. 

In Sec. g the clusters used in this work are described 
together with the basis sets and the methods of the cal- 
culations. We then consider in Sec. g first a (CuOs) 8 ~ ion 
and compare the results with predictions for a Cu 2+ ion 
in a single-electron approach. The results for the EFGs 
are given in Sec. [| and compared to previous theoreti- 
cal approaches and to experimental values. The magnetic 
hyperfine interactions are discussed in Sec. g| where spe- 



cial attention is given to study the transferred fields. The 
mechanism of spin transfer is analyzed in detail and the 
calculated hyperfine couplings are compared with values 
derived from experiments. Sec. |6| contains a summary and 
the conclusions. 



2 Clusters and theoretical methods 

The general idea of the cluster approach to electronic 
structure calculations of properties which depend predom- 
inantly on local structures is that the parameters that 
characterize a small cluster can be transferred to the solid 
and largely determine its behaviour. One should expect 
reliable results for local quantities like the EFG or the 
magnetic hyperfine fields. It would be much more ques- 
tionable to calculate quantities which are determined only 
by a large region of the solid. 

It would be desirable to have a cluster that contains 
as many atoms as possible, but there are two computa- 
tional limitations to the cluster size: the available com- 
puter resources and the convergence of the minimization 
procedure. The rest of the solid must be simulated by a 
shell of neighbouring pseudopotentials and, beyond that, 
point charges so that the cluster experiences the correct 
electrostatic potential. It is necessary, however, that the 
results obtained should be checked with respect to their 
dependence on the cluster size. 

An important feature of the cluster approach is that 
the atom to consider should be in a position of the cluster 
which respects as much as possible its natural surround- 
ing. This means that the cluster should retain the "local 
symmetry" around that atom. Since we are interested in 
the EFG of both the planar copper and the chain copper, 
we used clusters of two different symmetries according to 
which copper atom we were interested in. The alternative 
of choosing a single cluster to look at those two ions is 
currently prohibitive. However both these clusters include 
nearest neighbour atoms of the Cu02 planes, the latter 
being incorporated to reflect their proven importance. 

All clusters used assume the low temperature ortho- 
rhombic structure of YBa2Cri307, with lattice constants 
a = 3.827 A, b = 3.882 A and c = 11.682 A as given in 
Rcf. jy|. We adopt the standard labelling of the copper 
and oxygen ions as defined in Fig. |l|. Using the planar 
Cu(2) atom as a common centre, the crystallographic a- 
axis points towards 0(2), the 6-axis points towards 0(3), 
and the c-axis points towards 0(4). 

The smallest cluster used (see Fig. ||) consists of one 
Cu(2) in the centre surrounded by four oxygen atoms 
(0(2) and 0(3)) of the Cu0 2 plane and one 0(4) (which 
is called apical oxygen). Surrounding that central part, 
there are four pseudopotentials in the Cu02 plane with a 
valence of +2 simulating the core electrons of Cu(2) and 
one pseudopotential for the chain copper, Cu(l), with a 
valence of +1, twelve pseudopotentials for the Ba 2+ ions 
above the plane and twelve pseudopotentials for the Y 3+ 
beneath the plane. These pseudopotentials screen the elec- 
trons of the cluster from the naked positive point charges. 
The cluster is labelled CuOs/CusYi^Ba^. In addition, 
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these 35 atoms are surrounded by more than 6500 point 
charges. The positions of some charges at sites far from 
the cluster centre were adjusted in such a way that the 
correct Madelung potential in the central region of the 
cluster was reproduced. 

The largest cluster (see Fig. g) with a Cu(2) in the 
centre comprises in total 26 atoms all supporting a basis 
set; five Cu(2), sixteen planar and five apical oxygens. This 
cluster is surrounded by a shell of eight Cu(2), five Cu(l), 
twelve Y and twelve Ba pseudopotentials. That leads to 
the Cu502i/Cui3Yi 2 Bai2 cluster. 

To investigate the chain copper EFG it has been nec- 
essary to construct another cluster where a Cu(l) ion is in 
the centre. In detail the cluster contains 16 atoms, all sup- 
porting a basis set; one Cu(l), two Cu(2), two chain 0(1), 
two apical 0(4) and eight planar oxygens (four 0(2) and 
four 0(3)). There are also four Cu(l), ten Cu(2), eight Ba 
and eight Y bare pseudopotentials. Fig. [| shows the geom- 
etry of this cluster which is labelled Cu30i2/Cui4Y 8 Ba 8 . 

Table [lj contains a concise description of each cluster. 



0(4) 




0(3) Cu(2) °(2) 



Fig. 2. The Cus02i/Cui3Yi2Bai2 cluster centred at a Cu(2). 
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Fig. 1. The CuOs/CusY^Bai^ cluster centred at a planar 
copper, Cu(2). Ions with a full basis set are drawn darkly 
whereas ions with a pseudopotential are drawn lightly. 
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In view of the large range of cluster sizes used in this 
paper we continue our policy |17] of using the same ba- 
sis set for each cluster, so that direct comparison can be 
made, and as before use the 6-311G basis set [[19]. In any 
case, as earlier work has shown |l7J], there is no signifi- 
cant change if the quality of this basis set is improved. 
Pseudopotentials were used |2G] on Ba 2+ , Y 3+ , and those 
Cu 2+ and Cu + ions indicated in Figs, p], 0, and pi The total 
number of atoms, electrons, basis functions and primitive 
Gaussian functions for each cluster is collected in Table 0. 



Fig. 3. The Cu30i2/Cui4Y8Bag cluster centred at a Cu(l). 



In contrast to our previous work pffl on YBa2Cu30y 
we present here spin polarized calculations using the Gaus- 
sian 98 software package [ ji9| . 

We showed in an earlier work |L7| that the density 
functional (DF) method, a model which includes corre- 
lation effects, describes the covalent Cu-0 bonding bet- 
ter than Hartree-Fock (HF) methods. Although the DF 
operator, equivalent to the Hamilton operator in the HF 
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Table 1. Compilation of the used clusters with regard to their constituents. 



atoms with a full basis set 



atoms with pseudopotentials 



cluster centre 

Cu0 5 /Cu 5 Yi 2 Bai2 Cu(2) 

Cu 5 02i/Cu 13 Y 12 Bai2 Cu(2) 

Cu 3 Oi2/Cui4Y 8 Ba 8 Cu(l) 
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method, is not uniquely defined, the different exchange 
and correlation functional available give effectively equiv- 
alent results. In addition several forms of the gradient cor- 
rection of these functionals are available. So to maintain 
a consistent calculational scheme, following on from our 
work fll7|| on La2Cu04, we persevere with the exchange 
functional proposed by Becke pl],p2| coupled with the cor- 
relation functional of Lee, Yang, and Parr p3j (specified 
by the BLYP keyword in the Gaussian 98 program). For 
the final analysis a special program was written |24| which 
used the integrals and electronic wave functions from the 
Gaussian calculations as data. 



Table 2. Compilation of the used clusters: N: number of 
atoms, E: number of electrons, B: number of basis functions, 
P: number of primitive Gaussian functions. 



3 Embedded CiiCK ion 



ion see 



As a starting point, we consider the (CuOs) 
Fig. m) embedded in an appropriate lattice of pseudopo- 
tentials and point charges in the YBa2Cu307 crystal. In 
Table [| we collected results for the Mulliken charge pu , 
the components of the EFG V", the Mulliken spin den- 
sity p, the isotropic term a,i SO , and the diagonal elements 
of the dipolar hyperfine tensor aJL, at the copper site. A 
comparison with the values in the La2Cu04 system [[17| 
shows that the lacking apex oxygen in the YBa2Cri307 
compound has not a great influence on these values. 

We can compare these results for the hyperfine cou- 
pling constants to the predictions of a single electron pic- 
ture of a Cu 2+ ion. Bleaney et al. J25| assign to a single 
3d a ,2_ y 2 hole a dipolar coupling constant of 
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u dip ~ 

and a core polarization of 



To define the notation we specify that the EFG com- 
ponents, V m , are given here in atomic units, i.e. in ea^ 



L cp 



<r- 3 ) 



- K (r- 3 ) 



(3) 



(4) 



= - e a 



-3 



with a parameter k = 0.26 ± 0.06. 

In this paper we have preferred to use ai so instead 
of a cp since we are calculating the total Fermi contact 

is 



-Ha/a|. The quantity q u = V u /\e\ then interaction. This coukTbe dominated by a cp but 
corresponds to —9.7174 x 10 21 in units of V/m 2 . 

A general hyperfine interaction tensor, with compo- 



nents af3, is denoted by f a " and has units of ajj , which 
importantly is independent of the underlying nuclear 
structure. 

The hyperfine interaction energy at nucleus k is then 
given by 

kp a p = h k r) h ^ faf 3 ^ 

where k "f and 7 e are the gyromagnetic ratios for nucleus k 
and the electron respectively. 

Wc follow the NMR literature in replacing our general 
f a P by a a/3 , for the Cu on-site term, b a @ and c a ^ for the 
transferred hyperfine fields when it refers to Cu and O, 
respectively. 

We can identify at least three contributions to a gen- 
eral f af3 : 



however is the correct usage for 
the model system proposed. 



-3 



f a/3 
J tot 



JisoV 



ce0 



_|_ fO-P , fOif 

~ r J dip ' >' so 



a 6 



(2) 



where fi SO is the isotropic (Fermi contact) term, /^ is the 
traceless dipolar term, and f"@ the spin-orbit coupling. 
The difference between on-site and transferred contribu- 
tions will be discussed in Sec. pi 



If our dipolar hyperfine coupling a z ^ = —3.327 a 
is inserted into Eq. (J3|) the estimated value for (r~ 3 ) is 
5.82 a^ 3 which is close to the expected value. Assuming 
that the isotropic hyperfine coupling is dominated by core 
polarization we get a,i SO = a cp = —1.51 aTj from Eq. (Q) 
which is reasonable when compared to our properly cal- 
culated value in Table g| of a,i SO = —1.97 a^ 3 . We have 
used Eqs. (||) and (Q) to get these estimates and there- 
fore imposed the single electron picture of the Cu 2+ ion 
on our cluster results. We can question what our estimate 
for (r~ 3 ) of 5.82 aTj 3 really means since we are in a co- 
valently bound molecular situation. In the atomic picture 
the larger (r~ 3 ) is, the more compact is the associated 
atomic orbital (AO) and vice versa. In a molecular situa- 
tion we might expect a smaller (r -3 ) due to delocalisation 
of the orbital when involved in a molecular orbital (MO) 
but this is too simplistic for the following reason. 

The operator r -3 will sample only the density close to 
the copper and so should be similar to the atomic value if 
the electron was not shared with other atoms. This means 
that the (r~ 3 ) value will be proportional to the proba- 
bility of finding the electron in that orbital. In theoretical 
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Table 3. Compilation of results for the Cu(2) in the smallest clusters: a (CuOs) 8- ion in the YBa2Cu307 system and a 
(CuOe) 10- ion in the La2Cu04 system. The Mulliken charge and the Mulliken spin density attributed to the copper atom are 
denoted by pu and p respectively. The components of the EFG, V", are given in atomic units. The isotropic term, ai SO , and 



aerfine couplings, a^ ip , 


are in units of a s . 














PM 


yxx yyy 


yzz 
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O'iso 


XX 

^dip 


a yy 


^dip 


YBa 2 Cu 3 7 1.111 
La 2 Cu0 4 1.167 


-0.582 -0.485 
-0.706 -0.706 


1.067 
1.412 


0.633 
0.667 


-1.969 
-1.782 


1.651 
1.764 


1.676 
1.765 


-3.327 
-3.529 



calculations this is simply the square of the MO coefficient 
of the copper AO. In this context Eq. Q3J) is probably best 
rewritten as 



dip 



' — ^atom 



(r- 3 ) c 



(5) 



where n atom is the orbital occupancy. This equation is 
very useful if, and only if, (V -3 ) retains its value in 

a molecular situation. Otherwise Eq. (0) will have two 
unknowns but only one item of experimental data. 

We have extracted (r~ 3 ) , from our calculations 

\ / atom 

and find that the value (« 8 a^ 3 ) is significantly larger 
than the corresponding value given for the Cu 2+ ion. This 
is not entirely unexpected because it has long been known 
that AOs are more contracted in a molecular situation if 
we use AOs as basis orbitals in MO theory. It was particu- 
larly noticeable for hydrogen atoms and also for d-orbitals. 
This effect will be amplified if bonding reduces the elec- 
tron density in the AO, that is, makes the atom more 
positive. Although this has a passing interest it certainly 
raises a cautionary note here when using Eq. (||) or the 
modified form (H) . Since no useful conclusions result from 
our \i"~ 3 ) = 5.8 a^ derived from Eq. (0) using our theo- 
retical value for a z d z , it follows that no useful conclusion 
can be drawn from using any experimental value of a z d z 
in the same way. 



4 Electric field gradients 

To investigate the electronic structure of the Cu02 plane 
we used the cluster Cu502i/Cui3Yi 2 Bai 2 which com- 
prises five Cu 2+ ions. Calculations were performed for 
all three spin multiplicities m — 2, 4, and 6. The low- 
est energy is obtained for m = 4 probably reflecting 
that the spin distribution is reminiscent of an antifer- 
romagnetic system. The corresponding energy differences 
E m=6 - E m=4 and E m=2 - E m=4 are plotted in Fig. | 
versus the product po(pi + Pi + Pz + Pi) where po denotes 
the Mulliken spin density of the central copper ion and pi 
(i = 1, . . . , 4) that at the four other copper ions. It is seen 
that the calculated energies are related to the product of 
the Mulliken spin densities according to 



En 



E m=4 = jpo{pi + P2 + P3 + Pi) 



(0) 



with j = 283 meV. This empirical relationship and its 
connection with the antiferromagnetic exchange interac- 
tion will be further discussed elsewhere together with a 
comparison of the analogous results for other materials. 




Po^i 

Fig. 4. Dependence of the calculated energy differences 
-Em=6 — E m =4, and Em=2—E m =4, on the Mulliken spin densities. 



In Table |j the calculated values of the EFG compo- 
nents at the central copper atom are reported. The ex- 
perimental values for the quadrupole resonance frequen- 
cies |29| have been converted into EFGs using the value 
of Qp^Cu) = —0.211 b for the nuclear quadrupole mo- 
ment |2j . Also given in Table [| are previous theoretical 
values obtained from band structure calculations p2jO| 
and from a spin non-polarized cluster calculation"! 16|| . 
The present results are in good agreement with the ex- 
perimental data. Our values slightly depend on the spin 
multiplicity which reflects the influence of the finite size 
of the cluster since for m — 4 the absolute value of 
the (negative) Mulliken spin density at the central cop- 
per ion is smaller than that at the four other copper 
atoms. This will be further discussed in the next section. 
A comparison of the EFG values obtained in the small 
cluster Cu05/Cu5Yi 2 Bai2 (see Table Q) and the values 
at the four non-central copper ions in the large cluster 
(V zz = 1.41) shows good convergence with respect to 
cluster size. It should be emphasized, however, that the 
experimentally determined value for the asymmetry pa- 
rameter (jj = 0.01) cannot be accounted for by theoretical 
calculations which are based on atomic sites determined 
by crystallographic data. Although the orthorhombicity is 
small ((& — a)/b = 0.014) an asymmetry is expected which 
is appreciably larger than that observed experimentally. 

It is noted here that the EFG at the copper nucleus re- 
sults from relatively large and almost cancelling individual 
contributions, as is illustrated in Table ra for the example 
of Cu with spin multiplicity m — 4. This demonstrates 
that the theoretical determination of EFGs is quite deli- 
cate, as a further detailed analysis Eq] will show. There- 
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fore, it is necessary to describe the total electron density 
of the atom as accurately as possible, which require the 
basis sets to be sufficiently large. (In our experience, the 
6-311G basis sets Q have proved adequate as has been 
shown for La2Cu04 where the dependence of the EFG at 
Cu on the basis set was discussed in detail p7|.) 

A comparison of the theoretical results with experi- 
mental values shows that the new spin-polarized calcula- 
tions here are distinctly better than our previous spin- 
unpolarized calculations |10] and the FLAPW calcula- 
tions [§J||l]]. The FLAPW values are only about 50% of 
the experimental values although the contributions from 
the 3p and 3d orbitals (with the exception of the 3d x 2_„2 
orbital) are almost identical with our current calculations. 
The exception is crucial since according to our calculations 
the 3d a ,2_, y 2 orbital is significantly involved in covalent 
bonding with the 2p CT orbitals on the nearest neighbour 
planar oxygen atoms. The bonding to the apical oxygens 
is much more ionic. Clusters are basically molecular ions 
surrounded by an environment of atomic ions represented 
by pseudopotential functions or, more remotely, point 
charges. The wave functions are determined by a molecu- 
lar orbital approach which is ideally suited to the study of 
systems which are significantly covalently bonded. Aug- 
mented plane wave methods are more suited to systems 
which are mainly held together by ionic or metallic bonds. 
For this reason we believe that the cluster calculations are 
more appropriate than the FLAPW method for the study 
of local properties in these crystals. 

The values of the EFG components at the planar oxy- 
gen sites 0(2) and 0(3) are collected in Table g together 
with previous theoretical values and experimental data 
for resonance frequencies which have been converted into 
EFGs using the value Q( 17 0) = -0.02556 b. The agree- 
ment between theory and experiment is good. 

The cluster CusO^/Cu^YgBag was used to investi- 
gate the electronic structure of the apical oxygen 0(4) and 
the chain copper Cu(l) with spin multiplicities m = 3 and 
m — 5. The experimentally determined EFG components 
at Cu(l) exhibit the peculiar feature that the asymme- 
try parameter r\ is almost equal to 1 which is not due 
to a particular symmetry of the lattice site. Our cluster 
calculations with m — 5 fail badly to reproduce this asym- 
metry but those with m = 3 give r\ — 0.995 as shown in 
Table \Q For the latter, the spin densities reside mainly 
in the two planes with small values on the Cu(l) and the 
0(4). Nevertheless, the interplanar exchange coupling is 
of general interest. A careful study of this problem would, 
however, require much larger clusters (at least 11 copper 
ions) which are part of our plans for future investigations. 

For the apical oxygen, 0(4), we obtain V zz — —1.61 
and r\ = 0.22 also in reasonable agreement with experi- 
ments ID (| V zz |= 1.22 and rj = 0.32) in consideration of 
the small, but carefully chosen, cluster. 



5 Magnetic hyperfine couplings 



Table 5. Contributions of on-site AOs with spin up (a) and 
spin down (/3) to the EFG in the z-direction for the central 
copper and their sum calculated with spin multiplicity m = 
4. The remainder gives the contribution from nuclear point 
charges and AOs centred at other atomic sites. 
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5.1 Isotropic contributions 

In the previous section, we have shown that our wave 
functions can reasonably well predict the EFGs at vari- 
ous sites. Therefore we feel sufficiently confident to calcu- 
late other quantities from these wave functions, focusing 
in this section particularly on the magnetic hyperfine in- 
teractions in the Cu02 plane. 

The magnetic hyperfine interaction tensor at a certain 
nucleus can be decomposed into an isotropic part D and 
a traceless dipolar part T, the latter being discussed in 
the next subsection. The former is given by the difference 
between the spin-up and spin-down density at the nuclear 
site R 



D(R) = 



\ m m' 



(R)l s 



(7) 



where the sums extend over all occupied MOs. 
In a restricted open-shell calculation, the sum 



E 



in 



^(R) 



a 



2 would exactly cancel the partial sum Yl m ' 
2 counting over the number of spin-down MOs 



only and D(R) would be determined solely by the remain- 
ing (open-shell) MOs with spin up. In this case, D(R) 
corresponds to the Fermi contact term which is produced 
by AOs that are spherically symmetric around R. For un- 
restricted open-shell calculations, as they are discussed 
here, all MOs with spin up differ from those with down 
spin. If there are no singly occupied MOs with dominantly 
s-like AO contribution, as for example in the (CuOs) 8 ^ 
ion discussed above, the quantity D(R) is interpreted as 
core polarization. In the present context, we have both 
kinds of contributions. In a first approximation, the MOs 
with mostly 3d x 2_ y 2-character are singly occupied and 
yield a core polarization contribution to D(Cu). Neigh- 
bouring copper ions produce in addition transferred hy- 
perfine fields which produce a Fermi contact interaction 
from the 4s-orbitals and a much smaller contribution from 
the inner 3s-, 2s-, and ls-orbitals by spin polarization, in 
other words, a core polarization. We therefore prefer to 
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Table 4. Comparison of the EFG components in atomic units at the central Cu(2) in the cluster Cus02i/Cui3Yi2Bai2 
other theoretical approaches and with experiment. 



with 



multiplicity 


■trxx 


yvv 


yzz 


r\ 


this work 2 

4 
6 




-0.778 
-0.596 
-0.689 


-0.382 
-0.642 
-0.528 


1.160 
1.238 
1.217 


0.341 
0.037 
0.132 


other theoretical results a 

b 
c 




-0.841 
-0.31 
-0.30 


-0.823 
-0.27 
-0.29 


1.664 
0.58 
0.59 


0.01 
0.1 
0.02 


Experiment d 




-0.643 ± 0.002 


-0.628 ± 0.002 


1.271 ± 0.002 


0.01 


a Ref. [ 


16 


|, "Ref. [ 


L2 


|, c Ref. [ 


14 


|, d Ref. 


26 


• 



Table 6. Comparison of the EFG components in atomic units at the planar 0(2) and 0(3) in the cluster Cu502i/Cui3Yi2Bai2 
with other theoretical approaches and with experiment. 





0(2) 


0(3) 


multiplicity 
this work 2 

4 
6 


Trxx 

-1.347 
-1.341 
-1.101 


yvv 

0.772 
0.764 
0.684 


yzz 

0.575 
0.577 

0.417 


0.15 
0.14 
0.24 


\rxx 

0.718 
0.763 
0.662 


yvv 

-1.229 
-1.351 
-1.098 


yzz 

0.511 
0.588 
0.436 


0.17 
0.13 
0.21 


other" 

6 


-1.21 
-1.393 


0.72 
0.849 


0.49 
0.544 


0.2 

0.22 


0.72 
0.852 


-1.22 
-1.415 


0.50 
0.564 


0.2 
0.2 


experiment 


=F 1-08 


± 0.65 


± 0.42 


0.2 


± 0.65 


=F 1-05 


± 0.40 


0.2 


a Ref. 


12 


|, "Ref. [ 


14 


•' 


Ref. [ 


29 


1- 



Table 7. Comparison of the EFG components in atomic units at the chain copper in the cluster Cu30i2/Cui4Yi2Bai2 with 
other theoretical approaches and with experiment. 



this work 



multiplicity 
3 



■trxx 

0.601 



yvv 
-0.603 



yzz 

0.001 



other 



experiment 

Hi 



0.767 ± 0.003 



-0.773 ± 0.003 

a Rcf. |l2|, "Ref. [JL4|, c Rcf. ]l6[, d Ref. |ic|], e Ref. [J26[. 



0.995 



0.69 


-0.76 


0.07 


0.8 


0.565 


-0.631 


0.066 


0.79 


0.754 


-0.733 


0.022 


0.94 


0.404 


-1.153 


0.749 


0.3 



0.006 ± 0.001 0.99 



call D(R) just the isotropic hyperfine density at the nu- 
clear site R. 

The cluster calculations allow us to scrutinize the on- 
site and transferred hyperfine field in detail. We first fo- 
cus on the analysis of -D(Cu) using results from calcula- 
tions with the highest possible spin multiplicity. Thus in 
the cluster Cu502i/Cui3Yi 2 Bai2 with m — 6 we have 
Z?(Cu) = 0.13 a^ 3 at the central copper nucleus, but 
-D(Cu) = —1.33 a^ 3 at the four other copper nuclei. For 
the small cluster CuOs/CusY^Ba^, we got D(Cu) = 
— 1.97 &g . Furthermore, we have replaced copper ions 
in the cluster Cu502i/Cui3Yi 2 Bai2 by pseudopotentials 
without basis sets. This allows us to determine D(Gu) 
when there are two nearest neighbour (NN) copper ions. 
The values for D(Cu) obtained in this way are plotted in 
Fig. [?j as a function of the number of NN. 

It is seen, that the isotropic hyperfine density is pro- 
portional to the number of NN. In the same way it is 
found that D(0) at the planar oxygen site is about 0.6 a^ 3 
when there is one NN and about 1.2 a^ 3 for oxygens be- 



tween two copper ions. The sensitivity of the calculated 
spin densities can be assessed by an inspection of the two 
values obtained for D(Cu) with NN = 1 or NN = 2 in 
Fig. pi The slightly different values are due to the small 
orthorhombicity: the hyperfine density at the central cop- 
per nucleus in a cluster with three coppers aligned along 
the crystallographic direction a amounts to —0.99 a^ 3 but 
to —1.06 a^ 3 for direction b. 

In the same way we can analyze the results of calcula- 
tions with lower spin multiplicities. A careful analysis of 
all results shows that the isotropic hyperfine coupling can 
well be explained by the following ansatz: 



.D(Cuj) = a iso p{Cui) + ^ 



j£NN 



(8) 



Plotting D(Cva)/p(Cui) against J2jeNN P( Cu j) / P( Gn i) 
should then result in a straight line with slope /3j s0 and 
intersection a>i S0 . This is done in Fig. |6j. 
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Fig. 5. Isotropic hyperfine density (in units of a^ ) at the cen- 
tral copper nucleus plotted against number of nearest copper 
neighbours. 



13 

o 



13 
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-10 



Ep(Cu.)/p(Cu.) 



Fig. 6. D(Cui)/p(Cui) plotted against J2jeNN p( Cu j) I 
p(Cui). The different markers correspond to symmetrically in- 
equivalent copper sites in the cluster which differ by their num- 
ber of nearest neighbours. 



Similarly the isotropic hyperfine density at an oxygen 
site can be written as 



D(0) = llso ]T P {Cu 3 ). 
jeNN 



(9) 



Figure |7| displays the obtained values in a plot similar 
to the copper case. Again, the quality of the ansatz can 
be estimated by noting that the data points fit quite well 
to a straight line with slope ji S0 . From these two figures 
it seems quite clear that the influence of further distant 
atoms can be considered small. To make things more quan- 
titative we have assumed the isotropic hyperfine density 
at copper and oxygen to depend as well on next nearest 




E p(Cu.) 



J 

Fig. 7. Isotropic hyperfine density at the planar oxygen nuclei 
(in units of a^ 3 ) plotted against Y^jeNN P(^ u j)- The three dif- 
ferent markers indicate different spin multiplicities: o denotes 
m — 6, □ denotes m = 4, and A denotes m = 2. 



neighbours (NNN) and further distant copper ions accord- 
ing to 

D(Cui) = a iso p(Cui) + (3 iso ^ /°( Cu j)+ 

jeNN 

+ PL E -°( Cu ^ ( 10 ) 

j&NNN 

at a copper nucleus and 

D{0)= liao J2 P(C*i) + TL E <°( Cu ^+ 

jeNN jeNNN 

+< E p( Gu ^ ( n ) 

jGNNNN 

at an oxygen nucleus. A least squares fit to all data yields 
a iso = -2.68±0.02, f3 iso = 0.72±0.01, and (3' lso = -0.03± 
0.02 for the constants concerning the copper nucleus and 
7lso = 0.95 ± 0.01, | i i80 |< 0.025, and | 7 £ |< 0.03 for 
those concerning the oxygen nucleus (in units of a^ ) . 

It should be noted, that the marginal contribution of 
NNN ions to D(0) is of particular relevance for the in- 
terpretation of the 17 NMR relaxation rates. To demon- 
strate this small contribution we plot in Fig. the values 
of -D(O) as a function of NNN copper ions. The implica- 
tions of this result will be discussed further in Sec. p. 



5.2 Dipolar hyperfine coupling 

The dipolar hyperfine coupling results from a spatial av- 
erage of 1/r 3 with wave functions: 

T^(R) = £ (V4(r) | A^(r R) | ^(r)) - 
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Fig. 8. Isotropic hyperfine density (in units of a^ 3 ) at the pla- 
nar oxygen nuclei as a function of the number of next nearest 
neighbour copper ions. 



with 



-E(v6wi^(r-R)iv6(r) 



3x a xP - S af3 | r- | s 



A a/3 (r) = 



(12) 



(13) 



At variance with the isotropic hyperfine density which 
for copper exhibits a substantial transferred contribution, 
T"(Cu) is dominated by the on-site term. The transferred 
field from NN copper ions is small as can be seen in Fig. 
where T zz (Cu) is plotted as a function of the number of 
NN. 



-1 



O 



-3 



-4 




Fig. 9. The dipolar hyperfine field at the planar Cu(2) nucleus 
T zz (Cu) (in units of a^ 3 ) as a function of the number of near- 
est neighbour copper ions. Note that there is no substantial 
dependence of T zz (Cu) on NN. 



Again, all obtained values for the different clusters and 
multiplicities turned out to be proportional to the copper 
Mulliken spin densities and the same kinds of ansatz that 
worked for the isotropic case could be used. The results of 
this analysis can be written as 



T"(Cm) = ai oP (Cu i )+0, 



J dip 7 j 

jeNN 



and 



jeNN 



P(Cu,) (14) 



(15) 



The fitted parameters a di ■, /^ , and 7^ are collected in 



Table |. 



5.3 Discussion of spin density distribution 

The capability of the cluster approach to control the spin 
multiplicities enables us to get insight into the details of 
on-site and transferred spin densities. In Fig. [H] the differ- 
ence of spin-up and spin-down densities is shown along a 
bond axis for multiplicities m = 6 (upper panel) and m = 
4 (lower panel) calculated with the Cu502i/Cui3Yi 2 Bai 2 
cluster. The double humps at the coppers reflect the fact 
that the dominant part of the spin density distribution is 
related to the projections of the MOs onto the 3d a; 2_ a 2 
AOs. The small differences in the heights are due to the 
finite size of the cluster. The contact density at the two 
outer coppers is negative (they have only one NN) but pos- 
itive at the central copper (with four NN). At the oxygens, 
the spin density resides mainly on the 2p CT AO. Again, the 
tiny peaks at the positions of the nuclei reveal the contact 
densities which for the two inner oxygen atoms with two 
NN is twice as large as for the two outer oxygens with 
only one NN. The "antiferromagnetic" alignment (m = 4, 
lower panel of Fig. [l(]) also exhibits finite size effects. Since 
there are four copper ions with spin up and only one with 
spin down, there is an excess of positive spin density in the 
cluster which prohibits an exact cancellation of p^ and p^- 
at the bridging oxygens. The smaller spin density distribu- 
tion on the central copper atom in the antiferromagnetic 
case may arise from the removal of the restrictions due 
to the exclusion principle for the central spin compared 
to the m — 6 case. Note, however, that the electronic 
structure and spin density distributions for all the clus- 
ters are in perfect accordance if the analysis is performed 
with the individual Mulliken spin densities p(Cu) as has 
been demonstrated in the previous subsection (see Figs, o 
and|). 

An extrapolation from the finite clusters to an ex- 
tended system therefore just requires the knowledge of 
p(Cu). We assume that the value calculated for the central 
copper atom in the large cluster, p(Cu)= 0.68, is a reason- 
able estimate. The resulting hyperfine coupling constants 
(a-iso — 0.68ai SO and, similarly, for the other couplings) 
will be given and discussed in the next subsection. 
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Table 8. The values of the fitted parameters of our ansatz ( |l4| ) for the dipolar hyperfine field at Cu(2) and ( |l5| ) for that at 
0(2/3) (in units of a^ 3 ). Note that the dipolar field at the copper nucleus is nearly axially symmetric and a di refers to the 

crystallographic c-direction. For the planar oxygens, yV (7^) denotes the component in the CUO2 plane along (perpendicular 
to) the bond axis. 



J" 



A 



r 

dip 






idi'. 



idi; 



-5.16 ±0.01 0.136 ±0.004 0.47 ±0.01 



-0.26 ±0.01 



-0.21 ±0.01 




O Cu O Cu O Cu O 

Fig. 10. Spin density distribution (in &g ) along the crystallo- 
graphic axis a in the cluster Cus02i/Cui3Yi2Bai2. The upper 
(lower) panel displays the results from a calculation with spin 
multiplicity m=6 (m=4). 



To illustrate the mechanisms of spin transfer, the ra- 
dial dependence of the difference between spin-up and 
spin-down densities at the oxygen and the copper is shown 
in Figs, [ll] and [l2| for the clusters CuOs/CusY^Ba^ and 
Cu502i/Cui3Yi2Bai2. Owing to overlap and covalence ef- 
fects the Cu 2+ ion shares its spin density with the ligand 
planar oxygens whose spin direction is parallel to that of 
the local copper moment. The transferred spin is mainly 
on the O 2p CT orbital with a polarization of 8.2% and polar- 
izes the s-orbitals such that the nucleus sees positive spin 
density. This can be seen in Fig. O (left panel) where 
the radial dependence of D near the oxygen nucleus is 
shown separately for the projections onto the Is and 2s 
AOs. The contact densities are i?2s(0) = 0.737 a^ 3 and 
-Dis(O) = -0.154 a^ 3 . If the oxygen is between two NN 
copper ions (right panel) the contact densities are approx- 
imately doubled to D 2s (0) = 1.543 a B 3 and I>i s (0) = 



-0.304 a^ 3 . The densities | V>L(°) ? and 



1>L(P) | 2 are 



almost identical with values | tf>i s {0) 

I V , 2s(0) | 2 = 6.79 a^ 3 . Defining the polarizations /, 

the s-like orbitals by 



== 141 a B 3 and 
of 



Bn S (0) = y|i s (0)| 2 /„ 5 



(16) 



we get f 2s = 1.3% and f ls = -0.013% with NN = 1 but 
/ 2s = 2.7% and f ls = -0.026% for NN = 2. 
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Fig. 11. Radial dependence (along bond direction in units of 
&,b) of D (in a.g ) subdivided into contributions from MOs with 
mainly Is and 2s AO character at a planar oxygen nucleus for 
the small cluster CuOs/Cu5Yi2Bai2 with NN = 1 (left panel) 
and for the large cluster Cus02i/Cui3Yi2Bai2 with NN = 2 
(right panel). 



The isotropic spin density distributions near the Cu 
nucleus are shown in Fig. |l2|, for NN = (NN = 4) in the 
left (right) panel. The corresponding numerical values are 
summarized in Table O. The hyperfine fields transferred 
from the four NN ions change the polarization of the 4s- 
orbital from negative to positive. 



5.4 Hyperfine Couplings 

The results for the magnetic hyperfine couplings for the 
planar oxygen are represented in Table [lOJ. The spin-orbit 
contributions to the hyperfine couplings are expected to 
be small in the case of oxygen. We neglect them and as- 
sume that c" ot — Ci SO + c" ip . The values are very similar to 
the ones calculated for La2Cu04 (see Ref. |Q) with the 
exception that now Cf 0t is somewhat larger than c^ ot . 

Experimentally, these couplings cannot be measured 
directly but they are derived from a combination of various 
data. The spin shift measured with the external field in di- 
rection i is proportional to c" ot and the static uniform spin 
susceptibility. The spin-lattice relaxation rates are propor- 
tional to the sum of the squares of c^ t in the directions or- 
thogonal to i. Whereas the overall agreement between the 
theoretical and experimental values is satisfying, it should 
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Table 9. Values of the s-like AOs at the Cu (in a B 3 ), the isotropic hyperfine densities (in a s 3 ), and polarizations f ns for the 
cluster Cu0 5 /Cu 5 Yi2Bai2 (NN = 0) and the cluster Cu 5 2 i/Cu 13 Yi2Bai2 (NN = 4). 





NN = 


NN = 4 


n 
1 
2 


hMCu)|" 

7300 
725 


flns(Cu) 
-0.0586 
-3.542 


-9.59 x 10~ 5 
-5.83 x 10" 2 


|V 


«(Cu)r 

7300 
725 


£>» s (Cu) 
-0.0561 
-3.569 


fns[%] 

-9.18 x 10~ 5 
-5.88 x 10^ 2 


3 


107 


2.506 


0.280 




107 


2.530 


0.282 


4 


2.47 


-0.884 


-4.27 




2.36 


1.261 


6.38 



Table 10. Compilation of the magnetic hyperfine coupling constants in units of a B for the planar oxygen extrapolated to the 
infinite crystal. 





C-iso 


c 11 

'-tot 


c tot 


c tot 


this work 
Experiment" 


0.646 


0.966 
1.18 ± 0.08 


0.469 
0.69 ± 0.08 


0.503 
0.81 ± 0.08 


a Ref. [ 
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Fig. 12. Radial dependence (along bond direction in units of 
as) of D (in a B ) subdivided into contributions from MOs 
with mainly Is, 2s, 3s, and 4s AO character at a planar copper 
nucleus for the small cluster CuOs/Cu5Yi2Bai2 with NN = 
(left panel) and for the large cluster Cus02i/Cui3Yi2Bai2 
with NN = 4 (right panel). 



be noted that the calculated anisotropy, ciot/ c tot — %, is 
larger than the experimental value of about 1.7. 

When comparing the theoretical values with those de- 
rived from experiments (see Table 10) it is seen that the 
latter values are generally higher. The experimental anal- 
ysis relies on the approximation that the total spin sus- 
ceptibility Xs is isotropic and on the uncertain assumption 
that it is equally distributed over Cu(2) and Cu(l) sites. 
(In Ref. ||| the value Xs = 27.5 x 10 -5 emu/mole divided 
by 2.96 appropriate for the YBa2Cu306.96 sample was 
used.) Johnston and Cho fell reported anisotropic spin 



susceptibilities Xs (xj) of 37.4 (33.5) xlO 5 emu/mole. 
Using these values and the spin shifts from Ref. pQ] the 



hyperfine coupling constants would be c tot 
ci t = 0.58 a s 3 , and c c tot - 0.61 a s 3 . 

Yoshinari et al. pfl] also found that the hyperfine fields 
did not appreciably change between samples with an oxy- 
gen content of x = 6.6, 6.8, and 6.96. This is in sharp 
contrast with recent measurements |32] in the La2Cu04 
system where large differences between the hyperfine fields 
in undoped and doped (0.15% Sr) compounds have been 
observed. 

The calculated values for the isotropic and dipolar hy- 
perfine couplings for the planar copper are given in the 
left part of Table [ll]. We cannot yet calculate from first 
principles the spin-orbit contribution which for copper is 
not negligible. The total values for the on-site densities, 
a" ot = a iso + a^ + a" , however, depend crucially on a" OJ 
partly accidentally as will be shown below. We therefore 
choose to adopt the estimates given in Ref. |3^] which 
are based on the picture of a single missing electron in 
the copper 3d-shell developed by Bleaney et al. ]25| . From 
perturbation theory the hyperfine spin-orbit coupling den- 
sities are given by a phenomenological spin-orbit coupling 
parameter A, the energy differences AE between the var- 
ious 3d-orbitals and (r~ 3 )- Using the "standard" values 
for these quantities which are obtained from a variety of 
combined experimental and theoretical studies, one ob- 
tains oL = 2.409 a^ 3 and a^ = 0.427 a^ 3 . Following the 
arguments given by Monien, Pines, and Slichter |3J| who 
estimate the uncertainties in A and AE to about 20% and 
adding 10% for the possible change in (r~ 3 ), we include 

the 30% bounds in the 
Table [n]. 

In comparing our calculated values with previous esti- 
mates and values derived from experiment, we first note 
that we also get a transferred anisotropic hyperfine cou- 
,..|,;,.i, io ig% Q f ^.jjp va i uc f 5 iso _ Thg va i ue 

is considerably smaller than that of 



errors 



given for a\ ot and a^ ot 



pling b|L which is 18% of the value of b 



for b iso of 0.49 a B 3 
0.71 a^ 3 obtained in the La2Cu04 compound 0. This 
is due to the buckling of the planar oxygens in the CuC>2 
plane [p6| . Experiments give constraints on various com- 
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Table 11. The hyperfine couplings for the planar copper atom and comparison with literature values (in units of a B 3 ). The errors 
assigned to the total on-site hyperfine couplings originate from the uncertainties of the spin-orbit couplings as discussed in the 
text. (Since the dipolar hyperfine interaction tensor is traceless and nearly axially symmetric, aj[ ip 



-2 a Mp &nAb di P 



2 dip J 



disci 


"iso 


a L P b 


dip 


a" a 1 - 

U/ SO U'SO 


a" 

a tot 


a tot 


this work —1.82 


0.49 


-3.51 0.09 


2.41 0.43 


-2.93 ±0.72 


0.36 ±0.18 


other" -2.05 


0.74 


-3.67 




2.09 0.29 


-3.64 


-0.08 


b 


0.58 








-3.12 


0.40 


c 


0.71 








-3.01 


-0.18 


d 


0.69 








-2.75 


0.50 


e 


0.62 








-2.52 


0.59 


a Ref. [ 


|, "Ref. [ 


53 


], c Ref. [ 


34], d Ref. [ 


5|, e Ref. (35 


1- 



binations of the coupling parameters. Measurements of 
the 63 Cu NMR frequencies in the insulating antiferromag- 
netic compounds YBa2Cri306 and La2Cu04 give simi- 
lar frequencies of about 90 MHz. This then determines 



aiot 



4(fc iso + bj lp ) | = 



1.92 a R if an effective local- 



ized magnetic moment of 0.66 /j-b is assumed. Our values 
from Table [ll] yield 1.42 a^, 3 . A second constraint comes 
from the anisotropy ratio R of the nuclear spin-lattice re- 
laxation rates for fields perpendicular and parallel to the 
planes which is found to be about 3.7. Assuming antifer- 
romagnctically correlated spins R is given by 



R = 



Hb tso + b% 



lit - 4 (l 



i- 4- b 1 - 
iso T UfHp 



(17) 



Inserting our calculated hyperfine couplings into 
Eq. (|l7|) yields values for R which are far too large. An 
inspection of Eq. (|l7J), however, shows that R depends 
sensitively on a^ ot . In the nominator both terms add since 

a Lt < 0- I n the denominator, however, a cancellation takes 
place since a^ t > 0, and the relative values of the terms 
are important. Due to the accidental near cancellation of 



a iso and ajfc (a iso 



''dip 



0.068 aj/), the value of aj 



is crucial for a^ ot , as has been shown by Stoll et al. 



for 



La2Cu04 where, in addition, the spin-orbit contributions 
to the hyperfine field at the copper nucleus are strongly 
influenced by Sr 2+ dopants. Therefore the calculated cou- 
plings are inadequate in predicting the anisotropy of the 
Cu relaxation rates until more reliable values for a^ are 
available. 

A third experimental constraint comes from the fact 
that the copper magnetic shift measured for fields per- 
pendicular to the plane is independent of temperature and 
therefore assigned entirely to contributions from the chem- 
ical shift leaving a vanishing spin shift 63 K~s ■ This can hap- 
pen if a\ ot ± 4(6i SO ± &L) = 0. Our calculated values in 

Table |ll] give —0.61 a^ 3 and a value of a so — 3 a^ 3 would 
be required for a perfect cancellation. It should be empha- 
sized, however, that 63 K~s is found to be strictly temper- 
ature independent for most of the cuprates irrespective of 
the doping level. This is a remarkable experimental fact 
and an accidental cancellation of the various terms for sev- 



eral substances and doping levels is completely unexpected 
fronr our theoretical calculations of hyperfine parameters. 



6 Summary and Conclusions 

The local electronic structure of YBa2Cu307 has been in- 
vestigated by first-principles methods. Using clusters of 
various sizes, which attempt to incorporate the inrportant 
features of the Cu02-plane, the EFGs at the planar cop- 
per and oxygen sites were calculated with all-electron cal- 
culations and sufficiently large basis sets. The calculated 
values are in good agreement with the experiments. An 
inspection of Table pL where we examined in detail the 
various contributions to the EFG at the copper nucleus, 
demonstrates that the traditional interpretation of EFG 
values using Sternheimer shielding and antishielding fac- 
tors is highly questionable. 

The magnetic hyperfine properties in the planes have 
been determined very accurately by comparing the results 
fronr calculations with different spin multiplicities. It was 
shown that the transferred hyperfine fields at both Cu(2) 
and 0(2/3) are effectively due only to the nearest neigh- 
bour copper ions, the contributions from further distant 
neighbours being marginal. This finding has important 
consequences since it shows that the assumption of Zha, 
Barzykin, and Pines pj has no microscopic justification. 
The variant temperature dependencies of the Cu(2) and 
0(2/3) nuclear spin lattice relaxation rates can then only 
be explained by antiferromagnetic fluctuations which are 
commensurate. It can be argued that the neutron scatter- 
ing which exhibits incommensurate fluctuations is operat- 
ing on a much faster time scale than the NMR experiments 
or is probing periodic arrays of domain walls. 

The calculated values for the hyperfine coupling con- 
stants are in reasonable agreement with the values ex- 
pected from experimental data. As the discussion in Sec. |] 
demonstrated, the interpretation of the data is made 
difficult by the uncertainty of the contributions of the 
spin-orbit interaction. In this respect, improved quantum 
chemical calculations of the spin-orbit coupling would be 
desirable. 

The initial studies for this work were done in the 
diploma work of W.O. Dijkstra whose careful analysis is 
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